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Possibility of observing associated production of charged Higgs and W boson in the framework 
of MSSM at LHC is studied. Both leptonic and hadronic decays of W boson are studied while 
the charged Higgs boson is considered to decay to a r lepton and a neutrino. Therefore two search 
categories are defined based on the leptonic and hadronic final states, i.e. £ r+_E™ ISS and jj r+i5J" ss 
where £ = e or [i and j is a light jet from W decay. The discovery chance of the two categories 
is evaluated at an integrated luminosity of 300 fh~ x at LHC. It is shown that both leptonic and 
hadronic final states have the chance of discovery at high tan/3. Finally 5<r and 3<r contours are 
provided for both search categories. 

q . I. INTRODUCTION 

The Standard Model of particle physics has obtained a magnificent success through the experimental precision 
tests in the last years. Despite its success, the main part of the theory, i.e. the Higgs mechanism and the question 
of mass spectrum has remained unexplored and out of reach by the recent experiments. One of the main goals of 
upcoming experiments especially the Large Hadron Collider experiment (LHC) at CERN is to search for the Higgs 
boson as the predicted particle through the Higgs mechanism. While prediction of a single Higgs boson in the SM 
appears to be incomplete and reveals theoretical problems such as the Higgs boson mass divergence when including 
radiative corrections, theoretical models beyond the SM appear to be attractive candidate solutions to the SM Higgs 
O ,■ sector problems. The supersymmetric extensions to the SM are one of these models which not only provide elegant 
solutions to the Higgs boson mass divergence by introducing supersymmetric partners for SM particles but also 
^"'i provide approaches for other issues such as the gauge couplings unification Ref. [l(| . The family of two Higgs doublet 
models (2HDM) which introduce two Higgs doublets instead of a single scalar Higgs particle are of special interest 
'— among the supersymmetric models. The minimal supersymmetric extension to the SM, the so called MSSM, appears 
as a simple, reliable and minimal model belonging to the 2HDM family which predicts five physical Higgs bosons 
two of which are charged. While the neutral Higgs bosons of MSSM might be hard to distinguish from neutral Higgs 
1^ boson of the SM, the observation of a charged Higgs is a crucial signature of theories beyond the SM. This is an 
qq | enough reason that this particle has attracted special attention in the last years in different High Energy Physics 
f*"*- ■ experiments and will certainly be probed at the LHC experiment at CERN. 

The current results on the search for the charged Higgs include the low mass search by the LEP Higgs Working 
Group which excludes a charged Higgs with m/jj+} < 80 GeV in a direct search in Ref. Q while the combined result 
of MSSM Higgs boson searches excludes a light charged Higgs with rri(H+) < 125 GeV in the m/, — max scenario in 
Ref. [|[. The high tan/3 region has been excluded by the CDF collaboration quoted in Ref. Q. 
At LHC a large number of searches have been carried out in CMS and ATLAS collaborations. These searches divide 
the parameter space into two regions of light and heavy charged Higgs. The light charged Higgs (m^H+) < 175 GeV) 
■ is produced from the top pair production while the heavy charged Higgs (mtH+) > 175 GeV) is produced through 
^ ' 99 tbH~ and gb — s- tH~ with a combination procedure described in Refs. [5[ and p. 

In Ref. a light charged Higgs has been studied in CMS in the leptonic final state i.e. H + — >• t + v, W — > Iv, 
with I — e or /i. In Ref. [§[, a heavy charged Higgs has been analyzed in CMS in the hadronic final state, i.e. 
H + — > t + v, W — > jj. Both analyses predicted a possible discovery for a wide region of parameter space at an 
integrated luminosity of 30/& -1 . The H + — > tb decay channel has been analyzed in CMS in Ref. jl[ however it was 
shown to be suffering from the large hadronic background. 

The ATLAS collaboration results are reported in Ref. These analyses include three main final states for the 

light charged Higgs i.e. the hadronic decay of the r lepton with leptonic or hadronic decay of the W boson, and the 
leptonic decay of the r lepton with hadronic decay of the W boson. The heavy charged Higgs has also been analyzed 
by the ATLAS collaboration looking at H + — > t + v and H + — > tb. These analyses provide discovery and exclusion 
contrours corresponding to integrated luminosities of 1, 10 and 30 fb^ 1 on a significant fraction of the parameter 
space. There are also early data searches such as Ref. [ll[ which predict possible extensions to the current DO 
results in Tevatron with data as early as a time needed to collect 200p6 _1 integrated luminosity with ^/s = 10 TeV 
collisions. Although both H + —> t + v and H + — > tb decays are considered for heavy charged Higgs search, the latter 
does not lead to a promising discovery chance leaving the attention to be paid to the tauonic decay of the charged 
Higgs as the main decay channel throughout the parameter space. Therefore in the analysis presented in this paper 
only H + — > t + v is analyzed. 

The associated production of charged Higgs and a W boson is a complementary search channel to the current charged 
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Higgs analyses. This channel is attractive in the transition region where the on-shell ti production switches to the 
gg — > tbH~ and gb — > tH~ processes with the possibility of off-shell production of top quarks. Since the signal cross 
section decreases in the transition region around the top quark mass, the experimental selection of signal events is a 
challenge in this region. However since the H + W~ process can occure in both light and heavy charged Higgs area, 
it can help increasing the signal sensitivity especially in the transition region. 

This channel has been extensively studied in the literature from theoretical and phenomenological points of view. In 
Refs. [l2[ and [l3| the main production diagram was shown to be the tree level s-channel and i-channel bb — > H^W^ 
while the sub-dominant production process i.e. the one-loop diagram gg — > H^W^ which proceeds through the 
heavy quark triangular and box loops has a negligible contribution especially in the low charged Higgs mass region 
and high tan/3 values. 

In Refs. [l4{ and [l5[ quark and squark loops contributions to the 77 ± H /=F production process were evaluated and 
shown to contribute a little to the gg — > H^W^ 1 process keeping the bb — > H^W 11 still dominant. In Ref. [l6[ 
different choices of soft supersymmetry breaking parameters with or without squark mixing were studied and their 
effects on the gg — > H^W 11 process were evaluated. Again although these scenarios can have sizable effects on 
gg — > H^W^ , they were found to keep the bb — > H^W^ process as the dominant process. The supersymmetric 
electroweak corrections were also studied in Ref. [171 ] and they were shown to be negligible in the high tan/3 region 
where tan/3~ 30. In Ref. [l8| the 0(a s ) QCD contributions including virtual corrections as well as real gluon 
radiations were calculated in MS and OS schemes and the correction to the tree level process bb — > H^W T was 
estimated to be approximately —15% for the low mass charged Higgs and high tan/3. In Ref. [l9| supersymmetric 
QCD corrections were calculated and compared with corresponding SUSY-EW corrections (Ref. [13]) and 0(a s ) 
QCD corrections (Ref. 18]). As a result as an example at tan/3=40 it was shown that SUSY-QCD corrections are 
positive and can decrease 0(a s ) QCD and SUSY-EW corrections by ~ 50% for m^ H +^ — 200 GeV. 
The NLO QCD corrections to the cross section of this process were also calculated in Ref. [2(| and the results 
indicate that the K-factor (defined as (tnlo/&lo where <Jnlo(&lo) is the total cross section with NLO(LO) QCD 
corrections included) is close to unity for tan/3=40 and m^ H +^ = 200 GeV. Possible enhancement of the production 
cross section has also been studied in Ref. [2l[ in a general 2HDM and results show that there are regions in the 
parameter space where the 2HDM predicted cross section could be larger than that predicted by MSSM by two 
orders of magnitude. 

There have been several phenomenological analyses of this channel based on Monte Carlo simulation and event 
selection in the literature. In Ref. [22| the charged Higgs decay to a top and a bottom quark through H + — > tb was 
analysed and it turned out that there is no chance of signal extraction from the large hadronic background with this 
decay channel. Alternatively in Refs. 23 L 25] the H + — » t + v t decay channel was studied in the hadronic final state 
i.e. with W — > jj and the W + 2 jets process was taken as the main background in the analysis and concluded that 
in parts of the parameter space it is possible to have a signal significance exceeding 5a. 

The current analysis aims at three main purposes complementary to the previous analyses: first including a set of 
main background processes in the analysis i.e. ti, WW and JU+jets and estimating their contribution to the signal 
region, second extending the search to the leptonic final state which involves leptonic decay of W boson (i.e. W — > iv 
where I = e or //) and third extending the search to the light charged Higgs area where m H + < 175 GeV. The analysis 
starts with heavy charged Higgs which is studied in two categories of leptonic and hadronic final states. Then the 
light charged Higgs signal is analyzed and two production processes, i.e. H ± W T and ti — > H^W^bb are compared 
in both final states. A set of simulation packages is used for the analysis and a set of selection cuts is applied on 
the main distinguishing kinematic variables to increase the signal to background ratio. Finally the signal statistical 
significance is evaluated as a function of the charged Higgs mass and tan/3 and 5ct and 3<r contours are plotted. 



II. SIGNAL AND BACKGROUND PROCESSES 

The signal process which is analyzed in this work is bb — > H^W^ followed by the decay — > t^v t . The charged 
Higgs decay to tb is not studied in this work because it is expected to be hard to be distinguished from the large 
hadronic background although it starts to be the main decay channel with increasing charged Higgs mass in the 
region above the top quark mass. This argument follows the results quoted in Refs. § OH [22|- The charged Higgs 
decays to cs, c6, fj, + D and supersymmetric states are not studied in this work because they are negligible decays for 
heavy charged Higgs. For the light charged Higgs these decay channels grow but still remain at least two orders of 
magnitude smaller than the main decay channel i.e. H + — > tv Ref. [7[. The first two i.e. H + — > cs and H + — > cfe, 
suffer also from the large hadronic final state background. 

The reason for using only bb initiated process is that the dominant contribution comes from bb and the relative 
contribution of gg process becomes small and negligible (O(10~ 2 )) at high tan/3 which is the region of interest in 
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Signal, H^W* 


m (H ±) 


150 GeV 


175 GeV 


190 GeV 


200 GeV 


Cross Section 


0.15 pb 


0.12 pb 


0.103 pb 


0.095 pb 



TABLE I: Signal cross sections with tan/3=30. 





Background 


Process 


W+W~ 


tt 


W+jets 


Cross section 


115.5±0.4 pb 


878.7±0.5 pb 


187.1±0.1 nb 



TABLE II: Background cross sections used in the analysis. 



this study. This argument is in accord with theoretical studies of this process in Refs. |12[ and [13j. Depending on 
the W boson decay, two search categories are defined as the leptonic final state i.e. W — > Iv where I = e or \i and 
the hadronic final state i.e. — > jj where j is a light quark initiating a jet. Therefore the two signal processes 
under study are bb — > H ± W T — > T^v T lv and bb —> H ± W T — > T^v T jj. 

The background processes which are studied in this analysis are ti, WW and VF+jets. In the search for the leptonic 
final state, the tt process is forced to produce the following final states tt — > WbWb — > £TbbE™ lss or ljjbbE™ %ss . 
In the same search the following final state for WW process is considered: WW — > £tE™' iss or ljjE™ %ss . The 
VF+jets is used with W — > iv. All these processes are included to account for non-identified jets which escape the 
reconstruction and also the r jet fake rate which comes from light quark jets present in the event which pass the r 
identification algorithm. 

Subsequently the hadronic final state search is performed with forcing tt events to decay to the following final states 
tt — > WbWb — > rjjbbE'!p lss or TTbbE™ lss or jjjjbb while WW process is forced to decay to WW — > rjjE™ tss or 
ttE™ iss or jjjj. In the VF+jets simulation the W boson can decay as W — > tE™ iss or W — > jj. 



III. EVENT SIMULATION 



Signal events are simulated using PYBBWH code Ref. [26| which is linked to PYTHIA 6.4.21 Ref. [27|. The 
PYBBWH is currently the only Monte Carlo package for generating bb — » iJ ± M /T . This is the dominant production 
process in the high tan/3 region and provides a reasonable estimate of the signal rate. Therefore throughout the 
paper this process is used and declared as the signal process. As stated in the previous sections, although there can 
be negative contributions due to the higher order QCD or EW corrections including supersymmetric effects, these 
effects can be compensated by choosing a proper scenario of the 2HDM paremeters. 

The r leptons polarization and decays are controlled by the TAUOLA package (Refs. [28l - [30| ) which is linked to 
PYTHIA. Throughout the analysis r leptons are identified through their hadronic decays which results in a narrow 
r jet which will be described in detail later in the next sections. The analysis is based on parton showering and 
hadronization. The jet reconstruction is performed using PYTHIA built-in jet reconstruction tool, PYCELL, with a 
cone size of 0.5. Only jets within the pseudorapidity range of |?7| < 3.0 are reconstructed. The parton distribution 
function used in the analysis is MRST 2004 NNLO which is used by linking LHAPDF 5.8.1 Ref. [3l[ to the PYTHIA 
event generator. For the signal simulation the m^o — max scenario is used with the following parameters: = 200 
GeV, M- g = 800 GeV, n = 200 GeV and M SU sy = 1 TeV. 



IV. SIGNAL AND BACKGROUND CROSS SECTIONS 



The signal cross section is calculated at tree level by the joint PYTHIA+PYBBWH package using the MRST 2004 
NNLO PDF set. The background cross sections are calculated using MCFM 5.7 Ref. at next to leading order 
with the same PDF set as above. Table Q] lists the total cross section of the signal with different charged Higgs masses 
at tan/3=30 and Tab. [TT] shows the calculated cross sections of different background processes used in the analysis. 
The branching ratios used in the analysis are BR(t -> W ± b) = 0.98, BR(t — > H ± b) = 0.02 for mt H ±> = 150 GeV and 
tan/3=30, BR(iJ ± — » rv) = 0.987 which is almost constant for the mass range TOr#±) < 175 GeV. The W boson is 
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considered to decay according to the following branching ratios BR(W /± — > tv) = 0.106 and BR(H r± — > jj) = 0.68 
which are taken from Particle Data Group (PDG) at Ref. [33|. The branching ratios of charged Higgs decays are 
calculated using HDECAY 3.51 Ref. [Hj]. The branching ratio of top quark decay to charged Higgs is taken from 
PYTHIA. 



V. EVENT ANALYSIS AND SIGNAL SEARCH 

In this section the event analysis is described in detail. The search is divided into two regions of low mass and high 
mass charged Higgs with the separating point set to m/jj±\ ~ 175 GeV, which is the end point beyond which the 
production process it is pure SM process and no on-shell t — > H + b decay is allowed. 

The basic selection cuts are defined and used depending on the final state under study. In the leptonic final state 
a signal event contains a muon or an electron, from the W boson decay, with approximately the same kinematic 
properties, a r lepton from the charged Higgs decay and some amount of missing transverse energy originating from 
both W and r decays. Therefore a set of selection requirements is designed to identify these three physical objects in 
the event taking into account the kinematic differences between the signal and background to increase the signal to 
background ratio. Similarly a signal event in the hadronic final state contains two jets from the W boson decay, a r 
lepton from the charged Higgs decay and some missing transverse energy mainly from the r lepton decay. 
The kinematic thresholds applied on transverse momenta of muons and electrons as well as the missing transverse 
energy is proposed by comparing the distributions of these quantities between the signal and background samples. 
On the other hand the r identification algorithm used in the analysis tries to be as close to real LHC algorithms as 
possible. The leptonic decay of the r lepton is not analyzed due to the lower branching ratio of the leptonic decay 
and the soft lepton px spectrum from the r leptonic decay. 

In total four points in the parameter space are studied. One point is in the low mass region with m^H+) — 150 GeV. 
At this point two leptonic and hadronic final states are studied and compared with the main process which is the 
tt B^W^bb. The amount of excess one can observe by including the H^W^ signal is then estimated in the low 
mass region. 

The other three points i.e. rn^H+) = 175, 190 and 200 GeV lie in the high mass region and are studied each in the 
leptonic and hadronic final states separately. 

As will be seen in the next sections, the signal selection efficiency increases with increasing the charged Higgs mass 
due to the harder kinematics of physical objects in the event. However the cross section decreases more rapidly than 
the increasing rate of the selection efficiency, thus demanding higher tan/3 values to compensate the net decrease 
in the signal statistics when higher charged Higgs masses are studied. Therefore a 5<r and 3<r contour is plotted to 
present the accessible regions of the parameter space. Using these plots, a conclusion about higher charged Higgs 
masses can easily be made by extrapolating the contour. However the main conclusion is that the most sensitive part 
of the parameter space is the region studied in this analysis and other parts have less sensitivity and with increasing 
charged Higgs mass, higher tan/3 values are required to reach a 5cr discovery. 

A. Heavy Charged Higgs Search, Leptonic Final State 

In the following, a heavy charged Higgs with a mass of 175 GeV is analyzed in leptonic final state. This is a point 
in parameter space close to the defined border at which the on-shell tt —> H^W^bb process dies and the charged 
Higgs may only be produced by gg — > tbH~ + gb — > tH~ or ff ± V4 7=F . The aim of this work is to solely estimate the 
amount of excess one may observe over SM processes if a search is designed for H^W^ events. The set of selection 
cuts is described below for T[iE™ %ss final state, however a similar contribution from the electronic decay of the W 
boson i.e. TeE™ tss final state is expected. This is based on the assumption that in high energy events produced at 
LHC environment the mass difference between muons and electrons is negligible thus expecting the same kinematics 
between muonic and electronic final states and a reasonable lepton reconstruction efficiency achieved at LHC detectors 
ensures that no sizable lepton loss is expected. Discussion about the lepton isolation efficiencies and the performance 
of jet reconstruction algorithms in events containing electrons is beyond the scope of this work and needs a reasonable 
detector simulation. In the following event selection cuts are described for the final state containing muons, however, 
the signal significance is calculated assuming the same contribution from the final state containing electrons. An event 
must contain all the following objects to be selected: 
•One muon satisfying Eq. Q] 



p£ > 50 GeV, \rj\ < 2.5 
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FIG. 1: Transverse momentum distribution of muons in signal and background events in the T/i-EJ" ss final state . 



where r\ is pseudorapidity defined as 



V = ~ log 



tan — 
2 



in which 9 is the angle between the momentum and the beam axis. The adopted threshold in Eq. [T]has been inspired 
by Fig. m 

• One reconstructed jet satisfying Eq. O 

> 50 GeV, \r/\ < 2.5 . (2) 

This cut is also adopted in accord with Fig. [3] To reject the contribution of non-isolated muons which come mainly 
from heavy meson decays, the selected jet and muons with with pr > 20 GeV should be separated enough with the 
following requirement: 

Ai? (jct , M) > 0.4 (3) 



where AR is defined as AR = \J {Arf} 2 + (A</>) 2 . 

The r identification algorithm is then applied on the reconstructed jet in the event. A jet is accepted if it passes all 
the requirements of r-id. Although in a full simulation analysis there may be sophisticated requirements for r jets, 
here the basic requirements are applied. These requirements identify r jets with a reasonable purity. However the r 
fake rate (contribution of non-r jets which pass the r-id algorithm) may be higher than in a real analysis with full r-id 
cuts. This fact arises hopes for better results with less background contamination in the signal region. Verification of 
this issue is beyond the scope of this analysis which is providing prospects of the observability of this channel at the 
LHC and motivation of performing this analysis by LHC experiments, CMS and ATLAS. 

The r-id which is used in the analysis is similar to the one used by the CMS collaboration Ref. (35|. In this analysis 
all decays of r leptons are turned on but the search is only based on the hadronic decays which in total account for 
~ 65% of the total decay width of the r lepton. The r lepton in its hadronic decay produces predominantly one or 
three charged pions. Due to the low charged track multiplicity in the final state, the charged tracks (pions) in the 
r hadronic decay acquire relatively a higher transverse momentum compared to tracks of light quark jets. To study 
this effect, a jet-track matching cone of AR = 0.1 is considered around the jet axis. The hardest charged track in the 
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FIG. 2: Transverse energy distribution of reconstructed jets in signal and background events in the r/^i?™ lss final state. All jets 
with Et < 30 GeV are thrown away in the event as a general cut-off on the jet transverse energy. This is adopted usually in full 
simulation analyses to reject contribution of soft jets from pile-up events and also avoid reliability issues related to performance 
of the jet reconstruction algorithm at low transverse energies. 



matching cone is considered as the charged pion from the r lepton decay. Figure [3] shows distribution of the leading 
track transverse momentum in the matching cone around the jet. To exploit this feature a transverse momentum 
threshold is applied as the following: 

pleading track > 2Q ^ (4) 

Since r jets consist of few charged tracks, they are considered as isolated jets. To check if a r jet is isolated an 
isolation cone and a signal cone is defined respectively with a cone size of Ai? < 0.4 and AR < 0.07 around the 
leading track. The isolation is applied by requiring no charged track with px > 1 GeV to be in the isolation annulus 
defined as 0.07 < Ai? < 0.4. The low charged track multiplicity in the r jets also implies that the leading track in 
the jet cone carries a larger fraction of the r jet energy compared to quark jets in background events. This effect can 
be verified by plotting distribution of the leading track pr divided by the r jet energy as shown in Fig. HI In order 
to keep the signal statistics at a reasonable amount, a soft cut on this quantity is applied as the following: 

R = leading track pr / £7-7— jet > 0.2 . (5) 

The number of charged tracks in the r jet is also evaluated by a search in the signal cone and the following requirement 
is applied: 

Number of signal tracks = 1 or 3 . (6) 

Figure [5] shows distribution of the number of signal tracks in signal events before applying the cut. As seen from Fig. 
El r jets have undergone 1- or 3-prong decays predominantly. The azimuthal angle between the muons and r leptons 
is also investigated as shown in Fig. [51 As is observed from Fig. [51 background events tend to produce a harder 
back-to-back topology, however, in order to keep the signal statistics, no cut on this kinematic variable is applied. In 
the next step the r lepton charge is calculated as the sum of charges of tracks in the signal cone. Since muons and r 
leptons in signal events are produced with opposite charges, the following requirement is applied: 



Muon charge + r jet charge = 



(7) 
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FIG. 3: Transverse momentum distribution of the leading track in the r jet cone. 
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FIG. 4: Distribution of the leading track pr divided by the r jet energy. 
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FIG. 5: Number of tracks in the signal cone around the leading track of the r jet. 
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FIG. 6: The azimuthal angle between the muon and the r jet candidate in the T[j,E™ las final state. 
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FIG. 7: Missing transverse energy distribution in signal and background events in the r/^_B™ lss final state. 



•A reasonable amount of missing transverse energy should be left in the event. This is because of spin considerations 
and the fact that in iJ ± W /=F events, neutrinos tend to fly collinearly which is the result of charged Higgs boson being 
a spinless particle. The transverse missing energy distribution is studied as shown in Fig. [7] As seen from Fig [71 
the main backgrounds are all produced with a low E™ lss compared to signal events. The following requirement is 
therefore applied on E™ %ss : 

E^ ss > 50 GeV . (8) 

In order to have an estimation of the signal significance, all selection cuts are applied sequentially and their relative 
efficiencies with respect to the previous cut and the number of events which survive after each cut are listed in special 
tables. Table ITTT1 lists selection cuts efficiencies and the number of signal and background events remaining after each 
cut. The number of signal and background events remaining after all selection cuts (the last row of Tab. Mil multiplied 
by a factor of two to account for electrons) is used to calculate the signal significance according to Eq. [§] 

N s 

Signal Significance = —= (9) 
\/Nb 

Using Eq. [9] the signal significance with m/jj±) = 175 GeV and tan/3=30 is estimated to be ~ 0.4er with a data 
corresponding to 30 fb^ 1 integrated luminosity. At high luminosity run of LHC when 300/6 -1 data is collected the 
signal significance turns out to be 1.3a with m^u±) = 175 GeV and tan/3=30. The border of 5a contour (the point 
below which the signal significance is below 5a) starts from tan/3~ 60 and goes to higher tan/3values for heavier 
charged Higgs bosons. As a result a charged Higgs with a mass TO(h±) = 175 GeV is observable in the leptonic final 
state search if tan/3 > 60 and lower tan/3 values are out of 5a reach. 

It should be mentioned that using state of the art algorithms developed by LHC experiments, the fake rate is expected 
to be much smaller than what is observed in this analysis. Since lU+jets background is the main background and is 
a source of r fake rate, it is expected that this background would be under control in an LHC experiment analysis 
and lower tan/3 values would be in the 5a reach. 
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Process 




W + W~ 


tt 


W+jets 


Total cross section [pb] 


0.12 


115.5 


879 


187100 


Number of events at 30 fb~ 


369 


577577 


4394086 


5.9 xlO 8 


N Muons = 1 


172(46.7%) 


141662(24.5%) 


1.73e+06(39.3%) 


1.34e+07(2.3%) 


N Jets = 1 


76(44.2%) 


62159(43.9%) 


120121(6.9%) 


6.2c+06(46.6%) 


leading track pr > 20 GeV 


43(56.5%) 


19327(31.1%) 


35557(29.6%) 


1.5e+06(24.9%) 


Isolation 


35(81.7%) 


3676(19%) 


3247(9.1%) 


133281(8.6%) 


R > 0.2 


31(87.4%) 


2617(71.2%) 


2223(68.5%) 


86887(65.2%) 


1- or 3-prong decay 


30(98.4%) 


2177(83.2%) 


1916(86.2%) 


38389(44.2%) 


Opposite charge 


30(99.8%) 


2146(98.6%) 


1885(98.4%) 


35577(92.6742%) 


EJp isa > 50 GeV 


21(71.2%) 


776(36.1%) 


1112(59%) 


3737(10.5%) 


Total efficiency 


5.83% 


0.13% 


0.025% 


6.3xl0" 4 % 


Expected events at 30/6 -1 


21 


776 


1112 


3737 



TABLE III: Selection efficiencies and remaining number of signal and background events after each cut in the r/iZ?y" ss final 
state. The charged Higgs mass is set to 175 GeV and tan/3= 30. Numbers in parentheses are relative efficiencies in percent 
with respect to the previous cut. Branching ratios have been taken into account in transition from the second to third row. 

B. Heavy Charged Higgs Search, Hadronic Final State 

The heavy charged Higgs has been studied with this channel in hadronic final state in Ref . (23l - [25j j . In this analysis 
the intention is a closer look to this final state including main background samples and a r-id as close to what is in 
use by LHC experiments as possible. The event selection is described in the following. An event must satisfy all the 
following requirements to be selected: 

•There should be three jets satisfying Eq. [TO] The transverse energy threshold is set to the softer value of 30 GeV 
compared to the case of leptonic final state in order to avoid unwanted decrease of signal statistics when three jets 
are required. 

E^ > 30 GeV, \r)\ < 2.5 . (10) 
There has to be exactly three jets in the event and all pairs of jets are required to be well separated by Eq. [TTT 

AR UJ} > 0.4, Ai? (Tj) > 0.4 . (11) 

•The r-id algorithm is then applied on selected jets in the same way as is done in the leptonic final state search. 
Exactly one of jets must pass the r-id algorithm. Figures [5J [H [TU] show the distributions of related variables used in 
the r-id algorithm. 

•The other two jets remaining in the event are used for W mass reconstruction and their invariant mass is calculated 
as the W boson mass. The reconstructed W mass should fall within the W mass window defined in Eq. [T2"l 

\m(j tj) - 75.| < 20 GeV . (12) 

The adopted central value for the W invariant mass is a slightly lower than the current values. This is because the 
reconstructed jet energy may be different from the true value and the calculated invariant mass turns out to be slightly 
lower than the nominal value of 80 GeV. This is a known issue which is resolved by jet energy correction algorithms 
developed in LHC experiments. Figure [TT] shows the invariant mass of the two jets as the W boson mass candidate. 
•The missing transverse energy is also used to discriminate between the signal and background. Fig. [T^] shows the 
distribution of E™ zss in signal and background events. As seen from Fig. [T^the E™ zss distribution is much harder in 
signal events. Since the charged Higgs (H + ) is spinless, in its decay, it produces a left-handed v T and a left-handed 
r + . In the subsequent r + decay, the right-handed z7 T is kicked back to conserve the mother r + helicity. Therefore 
the two neutrinos preferably fly in the same direction producing a large amount of E™ lss in the signal event. The 
following requirement on E™ lss threshold is applied: 

E% lss > 50 GeV . (13) 

•A study of azimuthal angle between the r jet and E™ lss shows that background events tend to produce the neutrino 
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FIG. 8: Transverse momentum distribution of the leading track in the r jet cone. 
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FIG. 9: Number of tracks in the signal cone around the leading track of the r jet. 
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FIG. 10: Distribution of the leading track pr divided by the r jet energy. 
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FIG. 11: Invariant mass of the two jets as the candidate pair from W boson decay. 
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FIG. 12: Missing transverse energy distribution in signal and background events in the r jjE™ las final state search. 



back-to-back with respect to the r jet. The main source of E™ zss in background events is the neutrino from the W 
boson decay when W — > tv occurs. This decay prefers the r jet and the neutrino to fly back-to-back in the W rest 
frame. The change of the angle between the decay products when transferred to the lab frame is less than the case of 
charged Higgs boson and thus majority of background events tend to produce high values of A^ Ti£ m...). Since the 
charged Higgs is heavy (~ 2 x mw), when it is produced, it acquires a higher px compared to W boson (e.g. from 
W + W~ events) as shown in Fig. [TBI and thus gives a higher boost to the decay products leading to smaller angles 
in the lab frame. The difference is, of course, more obvious for heavier charged Higgs bosons. This feature results in 
lower A^( r £ ;mis S -) in signal events as shown in Fig. 1141 and is used by applying the requirement as in Eq. 1141 

Ac6 (TiBrss) < 2.5 . (14) 

With the above selection cuts, Tab. IIVI lists the selection efficiencies and the remaining number of events after each 
cut. Using the number of events as in Tab. |lV]the signal significance is estimated to be 0.5cr for rri( H ±- ) = 175 GeV 
and tanf3 = 30 at 30 fb^ 1 integrated luminosity. At high luminosity run of LHC, the signal significance could increase 
to 1.5(7. The tan/3 value for which the signal significance at high luminosity is 5tr is estimated to be tan/3 ~ 55. 



C. Light Charged Higgs Search, Leptonic Final State 



In the following, the contribution of H^W^ to the low mass region search in the two categories of leptonic and 
hadronic final states is estimated. In other words the additional contribution solely from _ff ± VF T to tt — > H ± W^ 1 bb 
which is the main signal in the low mass region is estimated. Having passed a common set of selection cuts, H^W^ 
signal should appear as an excess of events over what is observed from tt — > H^W^bb events. What is expected is 
that the final number of tt — > H^W^bb events should be dominant in the signal region due to the large total cross 
section of this process. One may imagine that a leptonic final state search for events with one jet (the r jet) in the 
final state may produce an iJ ± H A=F dominant sample. This is not the case because tt — > H^W^bb events have sizable 
contribution even to the one-jet bin due to their large cross section at the LHC (Fig. [T5]). To verify this fact the same 
set of selection cuts as what was used in the heavy charged Higgs leptonic final state search is applied on H ± W T 
and tt — > H ± W^bb events. The chosen parameters for the simulation are m^ H ±) = 150 GeV and tan(3 = 30. Table 



14 



Q. 

TJ 0.045 



~i — i 1 — i — I — i — i — i — i 1 — i — i — i — i — I — i — i — i — i — I — i — i — i — i — I — i — i — i — r 



p H± from hfW", m + =175 GeV 
pW* from W + W" 




150 200 250 300 

P T [GeV] 



FIG. 13: The transverse momentum distribution of charged Higgs with m^ H ±^ = 175 GeV in H events and W boson in 

W+W~. 



< 02 

^ 0.18 
■d 0.16 
^ 0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 




i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — r 



— Signal (H + W) 
■ ■■■ \ATVV" 
■ ■ ■ ■ tt 

W+jets 



miss 




Signal Final State: x jj E 
= 1 75 GeV 



■r i i i i ; 



".".. """"* "TLl I I t...]. 7 ' I I I I I I I I I I I I I I L_ 



0.5 



1.5 2 2.5 

A(Kt,E^ iss ) [rad] 



3 . 3.5 

miss\ 



FIG. 14: The azimuthal angle between the r jet and E™ " s in the r jjE™ tas final state search. 
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Process 




W + W~ 


it 


W+jets 


Total cross section [pb] 


0.12 


115.5 


879 


187100 


Number of events at 30 /b -1 


2363 


2154302 


16393214 


4.4xl0 9 


N Jets = 3 


652(27.6%) 


403831(18.7%) 


3.2e+06(19.5%) 


4e+07(0.9%) 


leading track pr > 20 GeV 


509(78%) 


195290(48.3%) 


2e+06(59.5%) 


1.8c+07(45.5%) 


Isolation 


318(62.5%) 


33839(17.3%) 


244890(12.9%) 


2.4e+06(13.2%) 


R > 0.2 


189(59.4%) 


12640(37.3%) 


95490(39%) 


752987(31.6%) 


1- or 3-prong decay 


185(97.8%) 


10104(80%) 


83556(87.5%) 


498813(66.2%) 


Exactly one r jet 


184(99.6%) 


10037(99.3%) 


83015(99.3%) 


497493(99.7%) 


W mass window 


98(53%) 


5543(55.2%) 


18803(22.6%) 


169107(34%) 


E^" ss > 50 GeV 


73(75.1%) 


2070(37.3%) 


12205(65%) 


22587(13.3%) 




57(78.2%) 


697(33.7%) 


5615(46%) 


7773(34.4%) 


Total efficiency 


2.43% 


0.032% 


0.034% 


0.000177% 


Expected events at 30/fo -1 


57 


697 


5615 


7773 



TABLE IV: Selection efficiencies and remaining number of signal and background events after each cut in the TjjE^p lsa final 
state search. The charged Higgs mass is set to 175 GeV and tan/3= 30. Numbers in parentheses are relative efficiencies in 
percent with respect to the previous cut. Branching ratios have been taken into account in transition from the second to third 
row. 
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FIG. 15: Jet Multiplicity in H ± W^ and tt — ► H ± W^bb events in the leptonic final state search with m^ H ±- ) = 150 GeV. 



[V] lists the selection efficiencies and the remaining number of events after each cut. Therefore as seen from Tab. [V] 
although the selection efficiency of the H^W^ process is roughly seven times larger they would appear as only 2% 
excess over ti — > H^W^bb events in the leptonic final state search with tan/?=30. 
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Process 






Total cross section [pb] 


0.15 


878.7 


Number of events at 30 


464 


108111 


N Muons = 1 


213(45.9%) 


42812(39.6%) 


N Jets = 1 


93(43.7%) 


11131(26%) 


leading track pr > 20 GeV 


50(54%) 


3284(29.5%) 


Isolation 


40(79.2%) 


1438(43.8%) 


R > 0.2 


35(87.7%) 


1214(84.4%) 


1- or 3-prong decay 


34(98.4%) 


1182(97.4%) 


Opposite charge 


23(99.9%) 


1158(98%) 


E^ las > 50 GeV 


15(66.4%) 


736(63.5%) 


Total efficiency 


4.9% 


0.67% 


Expected events at 30/6 -1 


15 


736 



TABLE V: Selection efficiencies and remaining number of H ± W zfl and it — > H ± W 11 bb events after each cut in the T(iE™* as 
final state. The charged Higgs mass is set to 150 GeV and tan/3= 30. Numbers in parentheses are relative efficiencies in percent 
with respect to the previous cut. Branching ratios have been taken into account in transition from the second to third row. 



Process 




it -»■ ^W^bb 


Total cross section [pb] 


0.15 


878.7 


Number of events at 30 /i> _1 


3020 


693544 


N Jets = 3 


800(26.5%) 


2e+05(30.3%) 


leading track pr > 20 GeV 


608(76%) 


1.4e+05(67.4%) 


Isolation 


367(60.4%) 


56938(40.2%) 


R > 0.2 


223(60.6%) 


31601(55.5%) 


1- or 3-prong decay 


217(97.4%) 


30463(96.4%) 


Exactly one r jet 


216(99.7%) 


30402(99.8%) 


W mass window 


116(53.7%) 


8664(28.5%) 


E^ las > 50 GeV 


81(69.4%) 


5441(62.8%) 




57(70.2%) 


3891(71.5%) 


Total efficiency 


1.87% 


0.56% 


Expected events at 30/6 -1 


57 


3891 



TABLE VI: Selection efficiencies and remaining number of H ± W zf: and it — » H ± W 11 bb events after each cut in the TjjE™ lsa 
final state search. The charged Higgs mass is set to 150 GeV and tan/3= 30. Numbers in parentheses are relative efficiencies in 
percent with respect to the previous cut. Branching ratios have been taken into account in transition from the second to third 
row. 



D. Light Charged Higgs Search, Hadronid Final State 

The contribution of H^W^ to the light charged Higgs search can also be estimated by applying the same selection 
cuts as was done in the heavy charged Higgs hadronic final state analysis. These cuts are applied on both H^W^ 
and it — > H^W^bb events for comparison. Table l*VTl shows the selection efficiencies and the remaining number of 
events after each cut. As seen from Tab. IVI1 _ff ± W /=F events appear as a ~ 1.5% excess over it — > H^W^bb events in 
the hadronic final state. Therefore in both categories of leptonic and hadronic searches, in the low mass region, the 
dominant production process is it — > H^W^bb. 

VI. 5a AND 3cr CONTOURS 

In order to obtain 5cr and 3tr contours, two more points in the parameter space are studied. They are rnm^) = 
190, 200 GeV with tan/3 = 30. Table IVIII shows the cross sections and efficiencies related to these points. It should 
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m (H ± ) 


Total cross section 


effjep. 


N lep. 


effhad. 


N had. 


190 GcV 


103 fb 


5.88% 


35 


2.74% 


52 


200 GeV 


95 fb 


6.04% 


30 


2.89% 


46 



TABLE VII: Cross sections and efficiencies of leptonic and hadronic final state searches for m^ H ±- ) = 190, 200 GeV and 
tan/3 = 30. eff; ep . (effhad. ) is the signal efficiency in leptonic (hadronic) final state search ( the leptonic final state includes both 
muons and electrons) and N l g P ' (Ng ad ') is the number of signal events which pass all the selection cuts at 30 fb -1 . 
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FIG. 16: The 5a discovery contour obtained with the leptonic final state search. The integrated luminosity is set to 300 fb' 1 . 

be noted that the branching ratio of charged Higgs decay to tv decreases with increasing charged Higgs mass. Using 
HDECAY, the following values are obtained and used in the analysis: BR(H ± — s- r^)=0.9 for m^ H ±- ) = 190 GeV and 
BR(H ± — > t^)=0.83 for rri( H ±^ = 200 GeV. It should be noted that although the signal significance stated by Eq. [9] 
grows like \[L where L is the integrated luminosity, in reality this assumption may not be the case especially when 
systematic uncertainties are introduced in the signal significance calculation. The estimation of such effects is beyond 
the scope of this work and needs a full simulation at the presence of detector effects. In this work discovery potential 
of the studied signal is presented. Figures [TBI and [T71 show the estimated 5tr discovery and 3a evidence contours if a 
leptonic final state search is carried out. Accordingly Figs. IT51 and [T9l show the ha discovery and 3a evidence contours 
with the hadronic final state search. As is seen the leptonic and hadronic final state searches turn out to provide 
almost the same search power for a charged Higgs boson with this channel and the available parameter space is almost 
the same for both final state searches. The excluded area is the extrapolation of CDF collaboration result in Ref. [3] 
to the heavy charged Higgs area and high tan/3 values. 

VII. CONCLUSIONS 

The associated production of charged Higgs and W boson were studied in the framework of MSSM. The signal 
observability in the leptonic final state search, i.e. 7? ± IU =F — » I tE™ %ss with I = e or fi and also the hadronic final 
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FIG. 17: The 3a evidence contour obtained with the leptonic final state search. The integrated luminosity is set to 300 fb 
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FIG. 18: The 5a discovery contour obtained with the hadronic final state search. The integrated luminosity is set to 300 fb 1 . 
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FIG. 19: The 3cr evidence contour obtained with the hadronic final state search. The integrated luminosity is set to 300 fb 



state search, i.e. H^W^ —> TjjE™ zss was investigated. The analysis is suitable for a high luminosity study at an 
integrated luminosity of 300 fb -1 . Results indicate that in the low mass region (m(#±) < 175 GeV) this signal only 
appears as a few percent excess (2% (1.5%) with mrjj±\ = 150 GeV and tan/3=30 in the leptonic (hadronic) final 
state) over what is observed from the dominant tt — > H ± W z ^bb process. In the high mass region at an integrated 
luminosity of 300 a charged Higgs boson with rri(H±) = 175 GeV, in the leptonic final state search, is in the 

5<7 discovery reach if tan(3 > 60. The 3<r contour starts from tan/3 > 46. The hadronic final state search leads to 5a 
signal for tan/3 > 55 with m^H ± ) = 175. The 3ct evidence in the hadronic final state search starts from tan/3 > 43. 
Although high tan/3 values are accessible with this channel, it can be used to extend the current Tevatron excluded 
region. Alternatively adding this signal to the current searches at LHC, especially in the high mass search, can provide 
broader regions in the parameter space than what is currently available and is also able to increase the signal statistics 
in the transition region. This motivates a full simulation at the presence of detector effects which can be done by 
CMS and ATLAS collaborations. 
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